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MODELING OF MACROSCOPIC STRESSES IN A DILUTE SUSPENSION OF 
SMALL WEAKLY INERTIAL PARTICLES 

ALEXANDER VIBE * 1 ’* AND NICOLE MARHEINEKE 1 


Abstract. In this paper we derive asymptotically the macroscopic bulk stress of a suspension of 
small inertial particles in an incompressible Newtonian fluid. We apply the general asymptotic 
framework to the special case of ellipsoidal particles and show the resulting modification due to 
inertia on the well-known particle-stresses based on the theory by Batchelor and Jeffery. 
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1. Introduction 

Suspensions of small, arbitrarily shaped particles in a Newtonian fluid are of great interest in 
physical, biological and engineering sciences, see among others 12 a Gang. Here one of the important 
questions, which often arises, is how the carrier flow is influenced by the suspended solids. This 
effect strongly depends on the particle size and mass, the characteristics of the flow but also on the 
scale of interest. In general, the initial configuration of the particles in the fluid domain involves 
some kind of stochastic nature, additionally the particles may interact with each other and also 
turbulent flow fluctuations may strongly alter the deterministic path of each particle. Thus looking 
on a cut-out of the domain on the microscale, where each particle has a significant size and that 
way a strong impact on the streamlines, the problem is dominated by stochastic effects. On the 
macroscale on the other hand, where the size of the particles is small compared to the dimensions 
of the flow, the suspension can be described as a homogeneous, non-Newtonian fluid and the task 
consists of modeling the corresponding bulk stresses. 

One step towards the derivation of the macroscopic suspension behavior is the study of the 
two-way coupled particle-fluid problem. Oberbeck described the disturbance flow generated by an 
ellipsoidal particle moving with a constant translational velocity through a stationary viscous fluid 
in na, Edwardes extended these results to the case of a rotating ellipsoid in [6j. In the case of a 
dilute suspension of small rigid spheres in a stationary Newtonian fluid, Einstein gave a derivation 
of the change in viscosity of the carrier fluid by first considering the motion of a single particle 
in a sufficient small region, such that the corresponding flow may be treated as linear, and then 
solving the stationary Stokes equations for the disturbance flow. The contribution of many particles 
followed by superposition of every single disturbance flow field [7. An analogous procedure was 
applied to ellipsoidal inertia-free particles by Jeffery [ID] . The model equation that describes the 
evolution of the principal axis of a prolate ellipsoid of revolution in a viscous fluid is known as 
Jeffery’s equation. Batchelor presented a general framework of modeling stresses in a suspension 
of rigid particles in an incompressible Navier-Stokes flow pQ. In the case of a dilute suspension, 
he applied his theory to ellipsoidal inertia-free particles using the results of Jeffery and derived 
an analytical term for the particle-stresses in the macroscopic description of the suspension. The 
study of non-dilute suspensions involves the modeling of the interactions between particles, see 
for example mill BUS]. Another general approach for suspension modeling was presented in [25] . 
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where Cauchy’s stress principle was applied to the suspension, reproducing classical results as well as 
showing additional effects induced by spatial non-uniformities of the dispersed phase. The influence 
of particle concentration was also considered in mm- Rigorous asymptotical homogenization 
strategies for suspensions were used for example in HIE]. The work by Junk & Illner m deals 
with the strict asymptotic derivation of Jeffery’s equation, using expansions in the small size ratio 
(ratio between the characteristic length scale associated to the particle and the one associated to the 
fluid). It yields a general strategy for constructing a correction to the undisturbed Navier-Stokes 
solution to account for the presence of a particle. 

In this paper we extend the asymptotic approach by Junk & Illner by taking into account inertia. 
We particularly classify three different types of inertia and use the asymptotic results for a single 
particle in the general framework of suspension modeling according to Batchelor. We study the case 
of a dilute suspension and consider particle-fluid interactions as well as gravitational forces. The 
regard of inertia requires the introduction of a special splitting for the forces associated with the bulk 
stresses. This way we deduce an extended influence of the solid phase on the macroscopic behavior 
of the suspension induced by small deviations of the particle movement from the streamlines of 
the carrier fluid (Theorem 1111) . Apart from a modification of the particle-stress tensor we obtain 
an additional effective particle force term. We demonstrate the inertia related effects for the well- 
studied example of a dilute suspension of (prolate) ellipsoidal particles fCorollarv fill). 

We organize this paper as follows: In Section [2] we set up the model for a single particle in 
a fluid flow with respect to different types of inertia. We give a brief overview of the asymptotic 
derivation including the essential definitions and assumptions. In Section [3] we sketch the basic ideas 
of Batchelor’s framework and use the asymptotic results to formulate corresponding consistent 
suspension models which account for weakly inertial particles. We illustrate the results for the 
special case of ellipsoidal particles in Section[J]and conclude with a summary in SectionjS] Technical 
details to the example are provided in the Appendix. 

Throughout this paper we use the following notation: 

Notation 1. Scalar values are typeset as ordinary characters a £ R. Vectors are indicated by small 
bold characters, v £ ]R", n > 1, their components are denoted by (v)j = Vi £ It, i = 1,..., n. In 
the following we will omit the ranges of the indices if they are obvious. Matrices are written as 
large bold letters, M £ R" xm with the components (M)jj = My, the identity matrix is denoted by I. 
Tensors of higher order are treated as linear mappings between the corresponding spaces and written 
as ordinary large characters, e.g. B : ]R 3 -A R 3x3 , or as the corresponding dyadic product. 

We use a tensor calculus notation, viz. v ■ w = Xu ViWi £ R for all vectors w, v, resp. (A ■ v)$ = 
V , ,l,,r, for A e R nxm , v £ R m , and (A ■ B) ik = V , A,for A <E R raxm , B e R mx£ . The 
dyadic product v <g) w £ R nxm of two vectors v £ R", w £ R m involves (v <g) w) • u = (w ■ u)v for all 
u £ R m , analogously for higher order tensors: Ui ® 112 ® ® u„. A two-dot product is denoted by 

(vCgiw) : (x®y) = (v-x)(w-y) /orv,x £ R n , w,y£ R m , analogously (u®v®w) : (x®y) = (v-x)(w-y)u 
and so on. For v, w £ R 3 the cross-product is given by (v x = Xu j e ijkViWj, with the Levi-Civita 
symbol 

{ 1, for (i,j,k) an even permutation of ( 1,2,3), 

— 1, for (i,j,k) an odd permutation of (1,2,3), 

0, else. 

We use the mapping B(y)ij = 'ffk^kjVk which assigns a vector to a skew-symmetric matrix such 
that B{y) ■ x = v x x for all x, v £ R 3 . 

As for the differential operators, the Jacobian is denoted by (<9 y f )ij = d Vj fi and the gradient by its 
transpose (V y f)y = d Vi fj. For scalar-valued functions these notations are certainly interchangeable. 
As usual, the nabla operator is formally treated as a vector and for example used for the divergence 
and rotation operators: div(u) = V ■ u. rot(u) = V x u. Higher order derivatives are noted as 9 y ... y u 
for a sufficiently smooth u. 


MACROSCOPIC STRESSES IN DILUTE SUSPENSION OF WEAKLY INERTIAL PARTICLES 


3 



Figure 2.1. Lagrangian description, bijective mapping between reference state 
and the actual time-dependent state. 


2. Mathematical model of a single particle in a fluid 

In this section we provide the ingredients for the suspension model. For this purpose, we consider 
a single small oriented particle suspended in a Newtonian fluid and study its influence on the carrier 
flow under the assumption that the particle is small compared with the dimensions of the flow. 
Proceeding from the three-dimensional interface problem we derive an asymptotic model for the 
particle in the fluid flow by help of expansions in the size parameter e (size ratio between particle 
and flow domain). The asymptotic analysis follows the procedure that was established for inertia- 
free particles in m- The novelty is the extension to inertial particles which we classify with respect 
to different types. 


2.1. Three-dimensional interface problem. 


Dimensional problem. A particle is a rigid body that we model as an open, bounded and time- 
dependent domain £{t) C 1R 3 with a smooth boundary d£(t) and time t £ Rg . The particle is 
characterized by its center of mass c : Rg" —> R 3 , c(f) = \£ (t)l ^ 1 ,[g( t \ xdx and an orthonormal right- 
handed director triad d; : Rg —> Sf, i = 1, 2,3, where Sf denotes the unit sphere in R 3 . Additionally 
we define a time-invariant reference state £ C R 3 with the center of mass lying in the origin and the 
directors identified with the Cartesian basis vectors {ei, 62 , 63 }. The function which maps every 
point of the reference state to the actual time-dependent state is completely presented by rigid body 
motions (translation and rotation), i.e. x : £ x Rg —» £(t), x(y ,t) = R(f) • y + c (t) (see Figure l2Tl) . 
with R € SO( 3) being the rotation matrix associated with the director triad, Rij = e, ■ dj. We 
assume that the small oriented particle is suspended in an incompressible Newtonian fluid. The 
regular domain f2 contains the fluid and the particle. The acting forces are due to particle-fluid 
interaction and gravity. The model of first principles for the particle and the flow consists for 
all t > 0 of the incompressible Navier-Stokes equations in the fluid domain S2 \ £(t), the no-slip 
condition on the interface as well as the kinematics and dynamics of the particle: 


p f (<9 t u(x, t) + (u(x, t) • V)u(x, t)) = V • S[u] r (x, t) + p f ge g , 

u(x, t) = u(t) X (x - c(t)) + v(f), 

^c(f)=v(t), ^R(t)=B(u(t))-R{t), 


(2.1a) 

(2.1b) 


V-u = 0, xgfl \ £(t) 
x e d£(t), 


(2.1c) 
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Pp\£ I^v(t) = J r S[u](x(y, £),*)• R(t) ■ n(y) ds(y) + p P \£\ge g , (2.Id) 

Pp\£\ ^ (J (t) ■ w(i)) = R(£) • J y x (R(i) T ■ S[u](x(y,f),f) • R(t) ■ n(y)) ds(y). (2.1e) 

Here, u and p denote the fluid velocity and pressure, S[u] = —pi + pj (Vu + Vu T ) the Newtonian 
stress tensor with pf being the fluid viscosity and p / the fluid density. The particle related quantities 
v and u) describe the linear and angular velocities of the rigid body. The inertia tensor is given by 
J(f) = R(t) • J ■ R (t) T with the time-invariant part J = \£\ _1 / e ||y|| 2 l — y<S>ydy and p p is the particle 
density. The gravitational acceleration is given by g with normalized direction e g . The system 
is completed by appropriate initial conditions for the flow and the particle as well as boundary 
conditions for the flow on dfl. 


Dimensionless formulation. We introduce two characteristic lengths associated with the fluid system 
x and with the particle y as well as a characteristic time t. We scale the velocities as u = v = x/t 
and the pressure with p = upf /x. For the sake of a simple notation, we keep the same symbols for 
all variables as in the dimensional formulation and get the dimensionless system: 

Re (d t u + (u ■ V)u) = V ■ S[u] T + ReFr -2 e g , V • u = 0, x€fi\£(t), (2.2a) 

u = u x (x — c) + v, x e 9f(f), (2.2b) 

^c = v, ^R = B(w).R, (2.2c) 

peRe |£|-^-v = / S [u] ■ R ■ nds + peReFr _2 |£|e„, (2.2d) 

at Jd£ 

pe 2 Re \£|— (J ■ lj) = R ■ / y x ("r 1 ■ S[u] R n) ds. (2.2e) 

at Jos ' ' 

The model m is characterized by four dimensionless parameters: the size ratio e = y/x, the 
density ratio p = p p /pf, the Reynolds number Re = uxpf /pf (ratio of inertial and viscous forces in 
the fluid) and the Froude number Fr = u/y/Wg (ratio of inertial and gravitational forces in the fluid). 
The Newtonian stress in dimensionless forixQ reads as S[u] = — pl+2E[u] with E[u] = 0.5(Vu+Vu T ). 
The size parameter enters the bijection between the reference and time-dependent state: 

x (y> t, e) = eR(£) • y + c(t), (2.3a) 

y (x,t,e) = e~ 1 R T {t) ■ (x-c (£)). (2.3b) 

The density ratio p is now the key parameter that allows us to define different inertial regimes. 


2.2. Tracer particles of different inertial type. As we will show our asymptotic approach 
results in particles which follow the streamlines of the surrounding fluid. Hence, we refer to them 
as tracer particles. To distinguish between different types of inertial particles, we introduce an 
e-dependent mass function a mQSS (e). With the particle density p p = p p (e) = cx rnass (e)m p /V p , the 
density ratio and the size parameter are then related according to 

m p Vf _ f ^ m p x 3 _ a mass (e) m 


_ Pp _ / \ 

P — — Qmassyt) T7 

Pf rnfVp 


3 (e) 


vl% 


7tt3 


\£\' 


(2.4) 


Here m pj V p and m/,b/ are the mass and the volume associated with the particle and the fluid, 
respectively, m = m p /mf. The freely selectable mass function has no physical meaning, but allows 
to balance the inertial terms in the particle momentum balance in different ways, yielding several 
models of inertial particles. In particular, we set a mass (e) = e k , k € IN, and define three inertial 
regimes with respect to the behavior of the density ratio in the limit e —> 0, i.e. heavy (k = 2), 
normal [k = 3) and light weighted (k > 3) tracer particles, see Table l2Tl For k = 3, the density 


*In the following S [.] always designates the Newtonian stress tensor of the corresponding solutions of the Navier- 
Stokes or the Stokes equations and E[.] denotes the associated deformation gradient tensor. 
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'- x mass c 

Name of type 

Behavior of density ratio 

k > 3 

light weighted tracer particles 

p —> 0 

k = 3 

normal tracer particles 

p = const 

k = 2 

heavy tracer particles 

p — OO 


Table 2.1. Classification of inertial types by means of the mass function, cf. m- 


ratio satisfies p = const leading to the inertia-free particle model that was investigated in m- The 
cases k < 1 are not covered by the asymptotic simplification as we will comment on in Remark [TJ 

2.3. Asymptotic analysis. In this subsection we generalize the asymptotic approach of m to the 
introduced inertial type classification. We particularly point out the differences in the asymptotic 
analysis that arise for the heavy tracer particles. 


Asymptotic expansion. We assume that under the considered inertia types the particle influences 
the carrier flow locally but not globally. This way we follow m and make the following assumptions. 

Assumption 2 (Particle-fluid interaction). 

1) The fluid is essentially undisturbed by the particle. 

2) The fluid motion induces a rotation of the particle of 0(1). 

3) The particle is far away from the boundary of the fluid domain dTl. 

For the quantities of the particle we take a regular expansion in powers of e, while the fields of the 
fluid flow are separated in a global and a local part: 

OO 

a = ^Va.; for a G {c, v,cj, R}, u = u 0 + R u /oc , p = Po + Pioc- 
i=0 

The local fields of the fluid quantities are expressed in the local coordinates of the particle reference 
state, 

OO OO 

u ioc(x,t,e) = y~yu; OCii (y(x,t, e),t), pi oc (x,t,e) = ^ e I ~ 1 p; oc ,i(y(x, t, e), t). 

i=1 i—l 

The scaling is chosen in such a way that the gradient of the pressure balances the Laplacian of 
the velocity. It is a consequence of the fact that the bijective mapping between the reference and 
time-dependent state is e-dependent, see (I2.3bl) . (This way the chain rule generates a factor of e -1 
every time u/ oc is being differentiated with respect to x or t.) 


Asymptotic solution. Our goal is to formulate a solution of the complete system (I2.2|) up to an error 
of O(e). Therefore, we consider only a finite number of asymptotic coefficients in the expansions. 
To address them we set 

2 1 

a e = ^e*a i foraGjc, v}, a e = ^^e*aj for a € {w, R}, (2.5a) 

2=0 

2 2 

u e = u 0 + R e ^Vu/ oc ,i, p e =p 0 +^V _1 £W,i- (2.5b) 

2=1 2=1 

Inserting the expansion coefficients into (12.31) yields the following form of the bijective mapping 
between the reference and the time-dependent state: 

x e (y, t, e) = eR e (^ e) • y + c e (t, e), y e (x, t, e) = e _1 (R e (t, e)) _1 • (x - c e (f, e)). 
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It is shown in ([TTj , Lemma 5) that (R e ) 1 = (R e ) T + 0(e 2 ) and || (R e ) 1 || = 0(1), thus y e is 
well-defined. The asymptotic coefficients of the mapping are particularly given by 
2 

x e (y,t,e) = ^Vxi(y,t), with x 0 = c 0 , xi = R 0 -y + Ci, x 2 = Ri-y + c 2 . 
i=0 

Moreover, we use Taylor expansions for functions of the form /(x(y, t, e)) in Co in the following, i.e. 

/(x) = /(c 0 ) + <9 x /(c 0 ) • (exi + e 2 x 2 ) + 0.5d xx /(c o ) : (exi + e 2 x 2 ) ® (exi + e 2 x 2 ) + (D(e 3 ) 

2 

= J2^f(c 0 ) + O(e 3 ), 

i=0 

with the differential operators 

D 0 = 1, Di = xi • V, D 2 = x 2 • V + 0.5xi ® Xi : V 2 . 

To emphasize the general structure of the asymptotic result in the subsequent Lemma [4] and to 
stress the relevance of certain terms for the suspension model in Section[3]we introduce abbreviations 
for some expressions. 

Abbreviation 3. Consider the following model-relevant functions for a ma ss(e) = e fc ; k >2 
hi = Rq • ( v i + B(u) o) ■ Rq • y — DiUo), 

h2 = Ro ‘ ( v 2 + (B(ujq) ■ Ri + B(u> i) ■ Rq) ■ y — D 2 Uo — Ri • u; OC) i), 

0, k > 3 
k 0 , k = 2’ 

{ 0, k > 4 
k 0 , k = 3, 
ki, k = 2 

gi = 0, k> 2, 

g 2 = -Ro • R 1 • gi - / y x (Rq • DiS[u 0 ] • Ro • n) ds + 

Jd£ v J 

where the derivatives of Uo and the Newtonian stresses are evaluated in Cq. Moreover, abbreviate 
the linear accelerations by 

-j - ( d I \ rj ■ d 

k 0 = to Re R 0 ■ I —v 0 - ^T 2 e s 1 , k: = TOReR 0 • —vi, 
and the angular accelerations by 


to 

= to Re Rq • 

d 

d t 

(Ro 

j • Ro u oj , 

tl 

= to Re Rj^ ■ 

d 

d t 

(Ro 

■ J • Rq 1 • lo\ + (r 0 • J ■ R'f + Ri • J • R^ • 


Lemma 4 (Asymptotic one-particle model). Let the following four requirements be fulfilled, then 
u e , p e , c e , \i e , u) e and R e defined in (1231) are a solution of the complete system (12.21) up to an order 
ofO(e) for ef 0. 

(RI) Let the global flow velocity Uo and pressure po be the solutions of the incompressible Navier- 
Stokes equations in Q x 1R + 

Re(<9 t u 0 + (u 0 • V x )u 0 ) = V x • S[u 0 ] T + ReFr -2 e g , V x • u 0 = 0 (2.6) 

with the boundary condition uo = u on dLl, and the initial condition uo(x,0) = u(x, 0) for 
x € LI \£(0), and u Q (x, 0) = w(0) x (x — c(0)) + v(0) for x G £(0). 


0, k > 3 
to, k = 2’ 



MACROSCOPIC STRESSES IN DILUTE SUSPENSION OF WEAKLY INERTIAL PARTICLES 


7 


(R2) Let the particle related coefficients satisfy the following set of conditions: Let 

— d(t) = Vj(i), i = 0,1,2 (2.7a) 

dr 

with the initial conditions Co(0) = c(0), Ci(0) = 0, i = 1,2 and let the zero-order velocity 
suffice the so-called ‘tracer condition’ 

\>o{t) = uo(c 0 (i),i). (2.7b) 

Lei additionally the matrices Ro, Ri solve the differential equations 

— Ro = B(u3q) ■ Ro, ^Ri = ' R o + ^( w o) • Ri (2.7c) 

with the initial conditions Ro(0) = R(0) and Ri(0) = 0. 

(R3) For every t > 0 let the local fields solve the stationary Stokes equations on the unbounded 
exterior of £ 


with the Dirichlet and integral conditions 


/ S[u* OCii ] • nds(y) = U, 
J d£ 

y x (S[u/ OCji ] • n) ds(y) = g,, 


Pioc,i fulfill the following decay properties: 


for some Ci, di, ff, * = 1,2 independent ofy. 

(R4) Let the mass function (12.41) be given as a mass {e) = e k with k > 2. 


— 0 

(2.8a) 

y e d£, 

(2.8b) 


(2.8c) 


(2.8d) 

iation [3J Additionally let 

U loc,i WTld 

, for y > fi> 0 

(2.8e) 


Proof. The proof results from a straight forward computation where the asymptotic coefficients 
are inserted into the full system m and Taylor expansions are applied for the global flow fields, 
whenever they are evaluated at the particle boundary (see |11] for normal tracer particles). The 
introduction of the different inertial types does not change the general mathematical structure of 
the problem but only affects the linear and angular momentum balances of the particle, i.e. the 
functions fq, g, in (R3). We sketch the steps of the proof for the sake of completeness. 

Inserting the asymptotic expansions for c, v, c o and R in (12.2cl) and applying (I2.7al) . (I2.7cl) shows 
that the particle kinematics (I2.2cl) is fulfilled up to an order of 0(e 2 ). 

The conservation of mass in (I2.2al) is fulfilled exactly since each local velocity field is assumed to 
be divergence-free. 

Considering the no-slip equation (12.2bl) . the use of the asymptotic expansions and the Taylor 
series of Uo in Co results in the following conditions that are exactly the Dirichlet conditions of 
uj oc ,i (I2.8bl) and the tracer condition (I2.7bl) . 

u 0 (c 0 (t),f) = v 0 , 

Diu 0 (c 0 (f), t) + Ro • uj oc> i(y, t) = vi + B(u 0 ) • R 0 • y, 

D 2 u 0 (c 0 (t), t) + Ri • u ;oCj i(y,t) + Ro • ih OCj2 (y,t) = v 2 + (B{u 0 ) • Ri + B(u i) • R 0 ) • y. 

Inserting the asymptotic expansion in the momentum equation of the fluid and applying the 
chain rule for the local fields yields 

Re(5 t u 0 + (u 0 ■ V x )u 0 ) - V x • S[u 0 ] T - ReFr _2 e s = R e ■ V y ■ (e _1 S[u ioCi i] T + S[u ioCi2 ] T ) + 0(e). 

Since Uo satisfies the Navier-Stokes equations and the local fields satisfy the stationary Stokes 
equations according to (RI) and (R3), the momentum balance holds up to an error of O(e). 
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As for the momentum balances of the particle, we insert the asymptotic expansion on both sides 
of the equations. On the right-hand side we use a Taylor series for the term S[u 0 ] in c 0 and keep 
the integrals for the local fields, this leads to 


m Re e fc 2 Rq • ( ^v 0 - ^e„ + e^-v 


At'" Fr 2 9 ' 'At 

= I S[u/ OCi i] • nds + e / (r,) • DiS[u 0 ] ■ R 0 + • Ri • S[u ioc ,i] + S[u; OCj2 ]j ■ n ds (2.9a) 


Ids Jd£ 

for the linear momentum and 

d 


mRee fc_1 R^ 


d t 


(Rq • J ■ Ro ■ u oj +(r 0 ■ J ■ Rq ■ w i + ^Ro ■ j ■ Rf + Ri ■ J ■ Rj) • 


!d£ 


y x S[u/ OCj i] • nd^ 


ids 


y x (Rq • DiS[u 0 ] • Rq + S[u /oCi2 ]) ■ nds + Rq ■ R : • J y x (S[u /oCil ] • n) 


(2.9b) 


for the angular momentum. These are exactly the integral conditions (12.8cl) and (I2.8cll) . 

The requirement (R4) is needed to avoid contradictory conditions as we explain in RemarkjT] □ 


In Lemma [H we assume that (R3) is fulfilled despite the fact that the problem seems to be 
overdetermined at first glance by the presence of Dirichlet as well as integral conditions on d£. This 
is however not true as the following lemma shows. 

Lemma 5 (Solvability conditions of Stokes problem). Let \M q , p q , q = 1, ..., 6 , be the solutions of the 
following six stationary Stokes problems related to the boundary conditions of the three elementary 
translations and rotations: 


VyS[w/ = 0, V y • Wq = 0, q = 1,.. -, 6, 
w/g = e q , w ?+3 = y x e g , q= 1,2,3, 

II W q || < Cg/||y||, HVyWglUbgll^dg/llyll 2 , 


yeR 3 \£, 
y e d£, 

llyll >fq> 0 


for constants c q , d q , f q . Set the following surface moments associated with w q , q = 1, 


( 2 . 10 ) 
6 to be 


Sg = / y x (S[wg] • n) ds(y), tg = / (S[w g ] • n) ds(y), 


IdS 


IdS 


Vg = / y ® (S[wg] • n) ds(y), W q = y ® y 8) (S[w g ] • n) ds(y). 
JdS JdS 

1) Let (R3) of Lemma^be fulfilled, then the following solvability conditions hold: 


Rq • ^0 
Ro ' V! , 


Ro W 
, R 0 • V 2 ; 


(' 


= ( Ro • ■ Ro ) : V 


):V, + ( 


Ro ■ <9 X Uq • Cl ) • tg 


eg-fi q = 1,2,3 
-eg_ 3 -g 1 <7 = 4,5 ,6 


( 2 . 11 a) 


(Rf ' o) + <9 x u^) • Ro — Ro" • (B(ujo) + 3 x U(f) ■ Ro • R^ ■ Ro + : Vg 

+ ( r.q ■ (Ri • Ro ■ (vi - <9 x u 0 ■ Ci) + 0.5[ci ® c x : V 2 ]u 0 + 5 x u 0 ■ c 2 )) ■ t g 
+ £ K Um {W q ) ktm 


e q ■ f 2 


k,£,i 


<7 = 1,2,3 
-eg_ 3 -g 2 9 = 4,5,6 


( 2 . 11 b) 


with (Ro) ij = e, • p , (L)ij = [p* ® c x : V 2 ](u 0 • p ) and K Um = 0.5[p fc <g> p f : V 2 ](u 0 • p m ). 
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2) Let ( 12 . 111 ) hold. Let {uioc,i-,Pioc,i) be the solution of the stationary Stokes problem ( 12 . 8 aD 
with the Dirichlet condition (I2.8bl) and with the decay property (I2.8el) . then the integral 
conditions (I2.8cl) and (I2.8dl) are fulfilled. 

3) Let (12.111) hold. Let (ui oc ,i,pioc,i) be the solution of the stationary Stokes problem (12.Sal) 
with the integral conditions (I2.8cl) and (I2.8dl) and with the decay property (12.8e[) . then the 
Dirichlet condition (I2.8bl) is fulfilled. 

4) The linear systems (12.111) are invertible. 


Proof. The key ingredient for this proof is the existence of the Green formula for solutions of the 
Stokes problem am Lemma 6 ), which connects the stresses and boundary conditions of two Stokes 
solutions and which needs the decay properties (I 2 . 8 cl) as well as those in ( 12 . 101 ) . 

For 1): The equations (I2.8cl) and (I2.8dl) are equivalent to 


&q ' fz 
e 9 • g, 


’dS 


S[ui OC}i ] ■ nds = w q ■ (S[u; oc ,j] • n) 
J dS 


ds, 


= e, • / B( y) ■ S[u/ OC; j] ■ nds = / (B( y) T • e q ) • (S[u ZoC)i ] • n) ds 


Ids 


IdS 


/ (y x e,) • (S[ui OC) j] ■ n) ds = - / w 9+3 • (S[u ioCii ] - n) ds 


( 2 . 12 a) 


( 2 . 12 b) 


J dS JdS 

for q = 1,2, 3. With the help of the Green formula the roles of w g and u; OCj i in the right-hand sides 
of (12.121) can be interchanged. Thus (12.121) is equivalent to 


e q ■ f, = / ^ • S[w q ] • nds, -e q ■ g ; = / h* • S[w g+3 ] • nds, q= 1,2,3. (2.13) 

JdS JdS 

We note that for i = 1,2 the terms h ; of Abbreviation [3] have the form 

= Rq ■ Vi + R„ ■ B{u>i_ i) • R 0 ■ y - = Rq • v; + B{ Rq ■ u?i_i) • y - iq, 

where we used M T • B(u>) ■ M = B{ M r • cu), which holds for any M € SO( 3) and any u> £ 1R 3 , and 
where iq = Rq ■ Diu 0 and r 2 = Rq • (Ri • uj OC) i + D 2 Uo — B(uj 0 ) ■ Ri • y). Inserting hi and sorting 
the resultant terms with respect to u>i_.i and v,; leads to 


(Ro •«*-!) 




■ S[w ? ] ■ nds(y) + 


* fz 

^q— 3 ' gz 


< 7 = 1 , 2,3 

q = 4,5,6’ 


(2.14) 


which directly results in ( 12 . 111 ) . 

For 2) we consider (12.111) which is equivalent to (12.141) and thus to (12.131) . since hi is the Dirichlet 
condition of Ui OCi i by assumption. We use the Green formula and get ( 12 . 121 ) . 

For 3) we analogously arrive at (12.131) . The assumption that the integral conditions are fulfilled 
is equivalent to (12.121) . Combining both results gives the equalities 



• (S[uzoc,i] • n) ds 


IdS 


hi ■ S[wq] ■ nds, 



w g+3 • (S[u; oc> i] ■ n) ds = 


h 

JdS 


h, ■ S[w g+3 ] 


n ds. 


Since Uz OCj i and w g fulfill the decay property the above equations imply u; OCj i = h.; on d£ by means 
of the Green formula. 

For 4): In ([TTj, Lemma 8 ) it is shown that the six vectors (t q ,s q ), q = 1,...,6 are linearly 
independent, thus the system ( 12 . 111 ) is invertible for all right-hand sides. □ 


Remark 6. 

1) The solvability conditions (12.111) together with the particle related equations (12.71) build a 
system of nonlinear ordinary differential equations (ODEs) for ci and R,_i, * = 1,2, where 
the ODE for the ith asymptotic coefficient depends on lower order coefficients, while the 
surface moments depend only on the geometry of the particle. The order of the corresponding 
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ODEs depends on the choice of the mass function a mass . For example in the case a m ass(e) = 
e, fi depends on dvi/dt and thus (12.1 laD becomes second order for Ci with an additional 
initial condition vi(0) = 0. However, under suitable regularity assumptions on uo, the 
corresponding right-hand sides of the ODE systems are at least continuous and thus at least 
local solutions exist by Theorem of Peano. 

2) Solving the solvability conditions of the Stokes problems (12.111) determine the Dirichlet con¬ 
ditions (I2.8bl) . However, it is not our aim to examine the existence, uniqueness and reg¬ 
ularity of the Stokes solutions given by (I2.8bl) , (12.8eD in general in this paper. Instead we 
will discuss some special problems in Section fjj where analytical solutions are available. 

3) Just like in [TT], Lemma [^] and\5\ imply a procedure for constructing a O(e)-solution of (12.21) : 
First the Navier-Stokes equations CUD are solved to get Uq and po, then the path Co of 
the particle is obtained from the tracer condition (I2.7bl) . The solution of the solvability 
conditions (12.11aD together with the corresponding differential equations (12.71) for i = 1 
provides the coefficients Ci, Vi and u>o, Rq that define the Dirichlet condition (I2.8bl) . The 
stationary Stokes problem (12.81) for i = 1 is solved. Finally, the last two steps are repeated 
to get u ; OCj 2 which completes the procedure. 

Remark 7 (Density ratio scaling). We want to point out the effect of the density ratio scaling 
introduced in m - The momentum balances of the particle CUD can be written as 

e fc_2 (k 0 + eki) = mi + em^, 

e fe-1 (€o + eli) = mf + em^, 

with k i, li of Abbreviation^ and m” +1 > m i+1 denoting the according integral terms appearing on the 
right-hand sides of (!2.9aD and (I2.9bl) . The terms ki,£i describe the linear and angular accelerations 
of the particle, while the integral terms connect them to the Dirichlet conditions of the corresponding 
local flow fields via the solvability conditions (|2. Ilf) . Different choices of the mass function lead to 
the following pattern for the balancing of the accelerations with the integral terms: 


£\k 

0 

1 

2 

3 

4 


-2 

k o = 0 

0 = 0 

0 = 0 

0 = 0 

0 = 0 


-1 

ki = 0 

k o = 0 

0 = 0 

0 = 0 

0 = 0 


0 

... = mi 

ki = m)’ 

pl-l 

E 

II 

o 

0 = mi 

0 = mi 

(2.15) 

1 

... = ml,’ 

... = mi 

E 

II 

ko = mi 

0 = mi 


2 

... = ... 

... = ... 

... = ... 

ki = ... 

k 0 = ... 


3 

£\k 

0 

1 

2 

3 

ki = ... 

4 


-1 

*o = 0 

0 = 0 

0 = 0 

0 = 0 

0 = 0 


0 

£i = m“ 

o 

II 

3 

^E 

0 = m“ 

0 = tn“ 

0 = m“ 


1 

... = m£ 

*i = m£ 

£ 0 = m£ 

0 = ini’ 

0 = rrij 

(2.16) 

2 

... = ... 

... = ... 

£i = ... 

£o = - 

0 = ... 


3 

... = ... 

... = ... 

... = ... 

£i = ... 

£o = ... 


4 

row shows 

the 0(e e )~correction of the momentum equations and each column corresponds 


to the choice of power of e of the mass function. 

For k < 1 the condition dvo/dt = Fr - ~ e g arises (see (12.151) ). In combination with the tracer 
condition (12.7bl) . which is independent of the mass function a maS s, this leads to the implication that 
duo(co(£),f)/dt = Fr -i e g in contradiction to the claim, Uo solving the Navier-Stokes equations in 
It, according to (Rl) of Lemma J This fact indicates that in this regime Assumption [H 1) does 
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not hold anymore, since the inertial effects of the particle are too strong and the perturbation of the 
surrounding fluid is of the same order as Uo . Thus, in consistence to the asymptotic, we restrict 
our classification of inertial types in Table\fffl ] to k > 2. 

2.4. Extension of Stokes solutions to the particle domain. The formulation of the suspension 
model requires the definition of velocities and stresses on the whole domain f 2 , including the interior 
of the particle. Hence, we need to introduce a meaningful extension of the quantities to the particle 
domain £. A possibility is to consider if as a part of the fluid domain and use the so-called singularity 
solutions which express the Stokes solutions by means of the Green’s dyadic (Oseen-Burgers tensor) 
m Such solutions have the disadvantage of not being bounded in the neighborhood of the origin, 
which is an important requirement in our modeling as we will see in Section [3j Alternatively, 
one could also think of a continuous extension of the Dirichlet conditions (I2.8bl) to £ and model 
the stresses as Newtonian. This leads to a discontinuity in the stresses at the particle boundary 
which contradicts with the desire for a smooth macroscopic stress on O. Therefore, we follow 
here the approach of [183 and treat the particle as a fluid with a rigidity constraint. Proceeding 
from a dimensional description (analogously to Section 12.11) . we obtain a consistent model for the 
dimensionless stresses and velocities in the particle domain by help of the presented asymptotic 
techniques. 

The fluid in the particle domain is characterized by the velocity z : £ (t) x R(J” —> R 3 of a material 
point and the non-Newtonian stress T : £(t) x R,} —> R 3x3 , i.e. 

p p (d t z + (z ■ V)z) = V ■ T t + p p ge g , Vz + Vz 2 = 0. x e £(t), (2.17a) 

z = u, T n = S [u] n, x € d£(t), (2.17b) 

for allt > 0 with the flow velocity u solving (I2.1al) and appropriate initial conditions in consistency 
with m■ The stress T acts as Lagrange multiplier to the rigidity constraint in (I2.17al) and can 
be expressed as symmetric gradient field of the three-dimensional unknown A : £(t) x R,} —»• R 3 , 

i.e. T = T [A] = VA + VA T [HjQ The velocity in the rigid body domain can be explicitly stated 

in terms of the particle quantities (center of mass c, linear and angular velocities v, cu of (12.11)1 as 
z = c)+v. Using this expression and the ODE for the director triad (12. lcl) simplifies (12.171) 

to a boundary value problem for T, resp. A: 

Pp ((^ R ) RJ ' ( x_c ) + ^ v ) = V ' T i X \ T + P P ( je g , x e £(t), 

T[A] • n = S[u] n, x £ d£(t). 

To transform the problem to the particle reference state we set T(x,f) = R(t) • T(y(x,f),f) • R 2 ( t ). 
In dimensionless form it is given by 

pR e f ^2 (e R ■ y + c ) - Fr -2 e g j = e -lR ■ V y ■ T , y £ £, (2.18a) 

T ■ n = R t ■ S[u](x(y, t, e), t) ■ R ■ n, y€d£. (2.18b) 

We model T as composition of the Newtonian stress of Uo and some disturbance, T(y,t, e) = 
R 2 (t) ■ S[uo](x(y, t, e), t) ■ R(t) + T; oc (y, t, e). As asymptotic expansion we take a regular power series 
in the size parameter e for T; oc and use a Taylor series in c 0 for S[u 0 ](x(y, t, e), t). With 

T = Tj OC; i + Rq ■ S[uo] ■ Ro + c {^ioc ,2 + Ro ■ -DiS[uo] • Ro + R 2 ■ S[uo] • Rq + R^ • S[uq] • Ri^ , 
we can formulate the asymptotic model for the stresses in Lemma [5] 


^The following considerations can equivalently be formulated for A, but since we are only interested in the stresses 
themselves, we do not use the inner structure of T explicitly It is however important to note that T have only three 
degrees of freedom such that the subsequent boundary value problems are well-posed. 
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Lemma 8 (Asymptotic model of stresses in particle domain). Let the requirements of Lemma [7] be 
fulfilled. Let T \ oc ^ = \7 y \i oc ,i + V y A J oci , * = 1,2 fee t/ie solutions of the following boundary value 
problems for t > 0 

Vy T L,i = f.. (2-19a) 

T; OC)i • n = S[u; OC)i ] • n, y € <9£, (2.19b) 

where the force terms are given by 

i [0, k > 3 

1 [l^ko, k = 2 5 

! 0, k > 4 

|£|-V fc = 3. 

|f|- 1 mReRS’-(^w 1 + ^ r R 0 . y ) ) fc = 2 

T/ien T , R e , c e sofee the system (12.1811 at least up to an error ofO(e). 


Proof. Similarly to the proof of Lemma [4j the approximation order follows from a straight forward 
computation. 

Inserting T on the left side of (12.18bl) . using the expansion described in the proof of Lemma [4] 
for S[u] on the right-hand side and taking (12.19bl) into account yields that (12.18bl) holds up to 0(e 2 ). 

When taking the divergence of T with respect to y, the terms of the form R, ■ S[uo] • R, vanish, 
since the Newtonian stresses are evaluated at Cq. With the asymptotic expansions of R and c, 
(12.18al) is fulfilled up to O(e) provided that (12.19al) holds. □ 


Remark 9 (Local velocity fields and forces in particle domain). 

1) By construction, the stress T induces a rigid body velocity z = e(dR/dt) ■ y + v in £, while 
its approximation T e = R 11 ■ T ■ R eT results in the approximation of z, namely 

d „ \ , / d „ 


z = v 0 


df 


Vl 


df 


V 2 


yg£. 


Since the functions Mi and R, do not depend on y, we can rewrite z e as 
z e = u 0 (x(y ,t,e),t) + eR e (hi(y,t) + eh 2 (y ,t)) + 0(e 3 ), 


using Abbreviationwhere hj(y, f) is the continuous extension of the Dirichlet conditions 
of Ui OC: i from Lemma ^ to £. This way we can extend the local disturbance velocity fields 
to £ by setting u; OC) j = h; OCi j. The so defined velocity is bounded by the assumptions of 
Lemma\4 } 

2) The force terms f, are not equal to fj of Abbreviation [21 but f F f^dy = holds true in 
consistency with (I2.19bl) . 


3. Kinetic model for a particle suspension 

In this section we present our main result, the asymptotical (macroscale) description for a particle 
suspension under consideration of inertial effects. For this purpose we set up mass and momentum 
balances by considering the particle suspension as a stochastic homogeneous fluid whose behavior 
is characterized by non-Newtonian stresses. To obtain a model for the qualitative behavior of this 
fluid, which is independent of the specific realization, we analyze the averaged equations. With 
the help of the ergodicity assumption we are able to derive an analytical expression for the bulk 
stresses. This procedure goes originally back to Batchelor in [I], the novelty of our work is the use 
of rigorous asymptotical results as model ingredients as well as the regard of inertial effects. These 
effects involve additional net volume forces and arise from the small deviations of the particles 
motions from the streamlines of the surrounding fluid. 
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3.1. General framework. Embedded into Batchelor’s framework [T] of suspension modeling, we 
introduce the necessary modifications to address the extra stresses induced by the particles’ inertia. 

We describe the suspension as a homogeneous fluid whose microscale properties vary stochasti¬ 
cally depending on the initial configuration of the positions and orientations of the immersed parti¬ 
cles. Given the probability space (W, A, p), let w . ,lo) : fix Rg" — > R 3 and g(.,., w) : fix Rjj" -a- R 
be velocity and pressure fields, u> £ W. The particle suspension is assumed to be described by (w, q) 
that is almost surely a solution of the incompressible flow problem 

Re (9 f w + (w ■ V x )w) = V x ■ Z r + ReFr -2 e CJ , V x ■ w = 0, for (x, t) £ fl x R + , 

supplemented with appropriate initial, boundary and E-related closure conditions. As in statistic 
turbulence modeling (see e.g. [ 23 ]) we decompose the random velocity field into a mean and a 
fluctuating part with expectation E[.] 

w(x, t, oj) = u(x, t) + u'(x, t, w), E[w] = u, 

analogously for the pressure q = p + p'. The task is to model the bulk stress Z or the inner force 
V x ■Z 2 , respectively. We particularly split the force into a divergence of surface stresses and a body 
force that is generated by the fluctuations of the translational velocity of the immersed particles 

V x ■ E T = V x ■ Z sT + b. (3.1) 

The viscous stresses are covered here by I s . The resulting averaged description of the suspension 
in fl x R+ is then given by 

Re (d t u + (u ■ V x )u) = V x • E[Z S ] T + E[b] - Re V x • E[u' ® u'] T + Re Fr" 2 e g , V x • u = 0 (3.2) 

(an analogue in turbulence modeling are the Reynolds-Averaged-Navier-Stokes equations). 
Assumption 10 (Suspension properties). 

(Al) The suspension consists of similar, independent and identically distributed particles £k(t). 

(A2) The suspension is dilute, i.e. the mean particle spacing is sufficiently big. 

(A3) Assumption [H (on particle-flow interactions) applies to every particle. 

(A4) No outer forces or moments act on the particles except of the gravitational force. 

(A5) The characteristic length scale of the fluid is much bigger than the length scale of the particle 
and the mean particle spacing. 

(A6) The fluctuations are generated by superposition of the local disturbance fields of the particles. 

(A7) The suspension is Newtonian up to an error of 0(e). 

(A8) The suspension is locally statistically homogeneous. 

(A9) The ergodicity hypothesis holds, i.e. for any involved random function f : flxR+xW —> R 71 , 
n > 1, there is almost surely an equality of the expectation and the integral of this function 
over a suitable space region V(x,f) C fl: 

E [/](M) = [ f{x,t,cj)dn(u}) = 1 [ /(£,*,w)d£ = E[/] v (x,f). 

Most of these assumptions were already used in [T]: With (Al) and (A2) particle interactions 
that arise if two particles are close together can be neglected. Assumptions (A8) and (A9) allow 
the switching from the abstract mean value to the analytically powerful volume average. To use the 
results of the asymptotic analysis we impose (A3)-(A7). In contrast to [I] we especially presuppose 
(A7) in order to consistently apply the asymptotical one-particle model of Section [2] Otherwise 
we would need to consider a non-Newtonian stress term of the surrounding fluid in the asymp¬ 
totics. Note that the introduced stress splitting 63 becomes essential for the proper handling of 
inertial particles when dealing with volume averages. Whereas the expectation E[.] is linear, the 
permutability of E[.] v and the divergence operator is generally not valid, see also Remark 1161 in 
Section [4] 
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3.2. Macroscopic stresses and forces. The core of the suspension model (13.21) are the macro¬ 
scopic stresses and forces that we deduce from the asymptotic one-particle model. We realize the 
fluctuation related quantities that are marked with the index ' by superposing the local disturbance 
fields -i oc of the particles (Assumption [TITl (A6)). 

Surface stresses. The derivation of the model for E s goes along pQ. It is based on the idea of applying 
the ergodicity hypothesis for E[E S ] and thus integrating the stresses over a suitable averaging 
volume V(x, £) which contains N(x,t) particles £k(t), k = 1,..., N. The underlying assumption on 
E[E s ] v (x, t) is 

E[Z s ] v = S[u] + _L r _p'| + v € u' + V € u' T d£+-L [ T'd£, (3.3) 

I v I Jv \uffc I v I Jue k 

where the stress fluctuations outside the particle are treated as Newtonian and inside the particles 
with respect to an appropriate extension. We particularly model the velocity and stress fluctuations 
by means of the local disturbance fields given in Lemma[4]for u' and in Lemma[8]for T\ With (A2) 

we may restrict to the influence of the local fields generated by the fcth particle on £k and neglect 

particle interactions. Using the two identities 

E / V« • d£ = E / {WiA + = / (V € ■ T' T )^ fc d^ + f T' k d£, 

j j &£ j j £g_ j £g_ j 

0 = E [d Xj u[) v = tL f %< d £ + T4 Z] / V « • (“S e j) d & 

\ V \ Jv\U£ e \r\ £ JSt 

and the divergence theorem on (13.31) results in 

E[Z*] V = S[u] - ^ [p' Id* + i4 E ( A T ' ■ ") ® € - u ' ® n - n ® ”' d ^) - /(V C ■ T' T ) ® ZdA 

\ V \ JV\U£ k \ V \k JS k J 

Hence, we obtain for the averaged surface stresses 

E[E S ] V = S[u] + S p , (3.4a) 

SP = M^(/ ( S [ Rfc - u L]-n)®^-R fc -uL®n-n®uf oc -R fcT d^) (3.4b) 

~ j £ (v c ■ (R fc • Tf oc ■ R fcT ) T )®^) 

where S p denotes the particle-induced stress tensor. The arising pressure term is incorporated here 
in the Newtonian stresses (as Lagrange multiplier to the incompressibility constraint). To express 
the integrals in S p (I3.4bl) with respect to the particle reference state we use the identities 

OO 

S[R ■ u; oc ] = E -e* _1 P*oc,J + e l <9 x (R ■ u; OCji ) + e l d x (R ■ u /oc>i ) T 
2=1 

OO OO 

= E £l_lR • (~P‘oc,i 1 + 9yU loc ,i + d y uf OCii ) R T = R E el_1 S[Uioc, l ] • R T , 
2=1 2=1 

J. oo 

V x • (r ■ T loc ■ R t ) = R E £l ~ 2v y • T L,i> 

2=1 



f(y(x,t,e))njds(x) 



e 2 E f 5 w/(y) d y = e2 [ /(y)( R • n )jMy) 

£ J & J dS 


and 
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for any scalar-valued smooth function /. This implies 


1 OO 


k 2—1 
i+2pfc 


'd£ 


s[uL,i] 


n ® (eR fe ■ y + c fe )di 


- e - ■ “K ■ I uf 4 ® n + n ® uf^dafr) ■ R" J - e l+i R B ■ / V y • T /ocl ® (eR K • y + c fc )dy. 
daf ' Je 

Using Lemma [SI the terms associated with the center of mass of each particle c k cancel each other 
out and, since f i is independent of y and J £ ydy = 0 by definition of center of mass in reference 
state, we get 


s ' = hE' :, r ‘ 


S[uf OCj i] ■ n ® (R K • y) - uf oc>1 ® n - n ® uf oc l ds(y) ■ R kT + 0(e 4 ). 


las 


Body forces. For the random body force we use a similar approach. We assume that the overall 
fluctuating force in the averaging volume V(x, t) is generated by contributions of each single particle 
in the corresponding particle domain 

b(x,t,w) = ^b fe (x,t,w)I ffc(t) (x), 

k 

where I denotes the indicator function, i.e. Ia(x) = 1 for x £ A and zero otherwise. Additionally, 
since the particles do not interact with each other, the random force generated by a single particle 
is modeled as the divergence of the local stresses in The average E[b] v (x,f) then follows as 

em v = p £ J v - ( R ‘ ■ T t • r‘ t ) t uns = ^ £ l - ( R * - T b ■ r ‘ t ) t 

= p E Rl - £ e ' +1 } v > ■ T ",. d « = jy £ Rt ■ £ t,+lf l+ 

k 2 —1 k 2 —1 

The last equality holds by Lemma [8] and Remark (9[ 2). 


Reynolds-kind stress tensor. The last step is the treatment of the Reynolds-kind stress term ap¬ 
pearing in & The key here is the fact that the local disturbance velocity is scaled with e and 
decreases fast enough at distance from the particle, while in its vicinity as well as in the particle 
domain it is bounded: By Lemma[4]the local fields fulfill ui oc ,i ~ (^(r -1 ) for r = ||y|| 0 and, since 

llyll = e_ 1 |l x — c||, it holds uioc ~ 0(e 2 ) for ||x — c|| 0. By (A2) and (A5) we find around every 

particle in V a ball B rk of radius r^, containing only the particle £k and it holds: 

E[u' ® u'] v = j^| J v u ' ® u ' d £ = |^| H Rk ‘ J v u tc ® u L d £ • R tJ 

4E Rfc 7 <c ® "LM ■ RkT + 0(e 4 ) 

l V U=l JB r k 

= lUT [ I u loc,i ® uf 1 dy + / uf 0 ! ® uf ol dy J ■ R ,,r + £>(e 4 ), 

fc=i \ J ^ k \£ ■’£ J 


where B? k denotes the ball centered at the origin with radius fk- The first integral in the last 
expression is bounded by Lemma|4]and the last by Remark[9] Altogether we get E[u'(g)u'] v ~ 0(e 4 ). 


3.3. Asymptotical suspension model. By means of the macroscopic stress and force descriptions 
we formulate an asymptotic suspension model to (13.21) that is consistent to the one-particle model 
and valid up to order 0(e 4 ) in the size parameter e. 
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Theorem 11 (Asymptotical model of a suspension with weakly inertial tracer particles). Let As- 
sumvtion \10\ be fulfilled, and let the local behavior of any particle in the suspension be determined by 
Lemma^4\and Lemma 0 Then the suspension is macroscopically described up to an error of 0(e 4 ) 
by 


Re (d t u + (u ■ V x )u) = V x ■ (S[u] + X p f 


Re Fr 2 e 0 


V x ■ u = 0 


m 


fl x R+, (3.5a) 


supplemented with appropriate initial and boundary conditions. The respective particle-induced 
stress and force are 


Z P (x,i) 


b p (x,t) 


1 

|V(x,i)| 

1 

|V(x,t)| 


N (x,f) 

E 

k=1 


e 3 R k 0 (t) 


N (x,i) 

E ( e 2 R o fc w 


k =1 


/ S[uf 1 ](y,t)-n®y 
Jae 

- hj(y,t) ® n - n ® h^(y,t)ds(y) • R£ T (t), 

•f?W + e 3 (Ro'W-f2W + R i'W-ffw)), 


(3.5b) 


(3.5c) 


with N(x,t ) being the number of particles in the averaging volume V(x,t) and with index k marking 
the asymptotic coefficients of the kth particle. 


The form of the particle-stress X p in (13.5bl) is known from the work by Batchelor [Tj, while 
our asymptotical approach for inertial particles (Section 12.21) give rise to an additional body force 
b p that is generated by small deviations of the particles’ center of mass from the streamlines of 
the surrounding fluid. Being 0(e 2 ), this extra-force dominates the particle contribution to the 
momentum of the fluid in the case of heavy tracer particles. One should further notice that the 
suspension behaves like a Newtonian fluid up to an error of 0(e 2 ) in consistency to the underlying 
Assumption [TUI (A7). 


4. Special case of a suspension with ellipsoidal particles 

The suspension of particles that have the same density as the surrounding viscous carrier fluid 
and the shape of prolate ellipsoids is often treated in literature; the structure of the particle-induced 
stresses is well-known in this case, see PQ and m■ In this section we illustrate the inertial particles’ 
effects by comparing our asymptotical suspension model of Theorem fill with those classical results. 
The determination of the respective analytical forms for X p and b p requires knowledge of the 
Newtonian stresses of the local disturbance velocity S[u; OCi i] at d £, an explicit form of the Dirichlet 
conditions for u/ OCi i and also of the differential equation for Ci, which is hidden in the solvability 
conditions in Lemma [5l Therefore, we first study the solvability conditions, deriving the ODE for 
Ci and an expression for u>o, also characterizing the Dirichlet conditions for ui oc ,i and then use the 
corresponding local stresses to evaluate the integrals in (13.51) . The necessary calculations are quite 
technical and lengthy as they are strongly related to the geometrical properties of the ellipsoids. 
For the sake of completeness they are summerized in the appendix. 

4.1. General properties. We start with introducing the quantities that characterize the geometry 
of arbitrary ellipsoids and discuss then the implications arising for the asymptotical framework. 

An ellipsoid is given by £ = D Bi, where B\ denotes the unit ball in 1R 3 and D = diag(di, d^, ^ 3 ) 
is the diagonal matrix with the lengths of the semi axes di > 0, i = 1,2,3. The surface moments 
appearing in the solvability conditions (12.111) in Lemma [5] can be provided as 

s <7 = 0. s g+ 3 = Cq+ 3 |^|(D 2 — tr(D 2 )l)e 9 , t q = 2,C, q \£\e q , t g+3 = 0. 

V g = 0, M q+3 = C q+3 \£\D 2 -B(e q ), W q = ( q \£\D 2 ® e q , W q+3 = 0, 

with the geometry dependent constants ( q , q = 1 , 2 ,3. Appendix lA.2l is dedicated to their derivation, 
the explicit expressions for f q are stated at its end. From (12.111) we see that the structure of the 
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surface moments implies a decoupling of the conditions for the linear and angular velocities, hence 
we get 


Ro ' v i = R o ' dxU 0 • ci + (3|£| diag(Ci,C2,C3)) 1 • fi, 

Ro • w 0 = 0.5Rq • V x x u 0 + |f| _1 (trD 2 l - D 2 ) -1 • diag«4,Cs, Ce) -1 • gi 


+ (tr D 2 I 


D 2 ) 1 • diag (d\ —d\,d\ —d\,d\ 


dl) 


((Ro • E[u 0 ] ■ r o) 32 \ 
(Ro • E[uo] ■ Ro)l3 

V R o • E[u 0 ] ■ Ro) 21 / 


(4.6a) 

(4.6b) 


with E[uo] = 0.5(V x u o + V x Uq ). In our set up, g : = 0 always holds (cf. Abbreviation [3|. Conse¬ 
quently, the Dirichlet condition (I2.8bl> is presented by 

hi = Rp (vi - <9 x u 0 • cr) + Rp • (B(u 0 ) - <9 x u 0 ) R 0 • y = h“ ns ‘ + Ai y. 


where the matrix-valued function Ai = A] (t) is independent of y and the time-dependent vector 
hj 0 "' st (t) becomes 

^const _ ( 3 |£.| di a g(C g , q = 1,2, 3)) _1 • f l5 

because Vi — 9 x up • Ci is exclusively determined by the source term in (14.6aj) . The last ingredient 
for the suspension model are the Newtonian stresses of uj OCi i at d£. Since the Stokes problems in 
Lemma[4]are linear, the local velocity field u; OCi i is given by the linear combination of the Oberbeck 
uioc,i,Ob and Jeffery solutions u/ OC)lj j e (see Lemma lA.21 stated in Appendix IA.il) . This yields the 
following stress terms that can be computed by means of the techniques provided in Appendix IA.21 


> as 


S[u ;oCjli oh] • nds = f 


i) 


'as 


S[u/ oc>lj j e ] • nds = 0, 


'as 


S[u; OCi i j0 b] • n ® yds = 0, 


<g> n + n <g> h™ nst ds = 0, 


(4.7a) 

(4.7b) 


/as 


S[u; OCj i iJe ] • n <g> y - Ai • y <g> n - n <g) Ai - yds = \£| -M- 4(M : diag(a°, a %, a^))! , (4.7c) 


J as 

where 8 = did^ds- For the definition of the geometry dependent scalars a° and matrix M we refer 
to Appendix lA.il 

4.2. Suspension model of weakly inertial ellipsoidal particles. We combine the results from 
the one-particle asymptotics and deduce the macroscopic suspension description for arbitrarily 
shaped tracer ellipsoids (TE) that we even specify for prolate ellipsoids in Corollary H4l 

Arbitrarily shaped ellipsoids. The particle-induced stress tensor for arbitrarily shaped ellipsoids 
becomes by means of (13.5bl) and (14.71) 


57 = e3 |E( i R o • M " • R o fcT - 4 (M fc : diag(a;, «§,«§))■) . 


N 


|V| 


k= 1 


(4.8) 


The force term b ?J (I3.5cl) depends via the functions ff (Abbreviation [3]) on the choice of the mass 
function (12.41) . With respect to the three different inertial types of tracer ellipsoids we obtain 

bP = e2 p^E( mRe -+ e (~|£| v • S[uf + m R e 


= 6 




^WS|uf + £ V(J|ReE 


d t 


Vl 


- v-sfuV 


0(e 4 ), 


heavy TE, 




















18 


A. VIBE AND N. MARHEINEKE 


b " = ^ E (-I f l V ' S M T + ) 

= e 3 (f> — 1^ V • S[u] T + 0(e 5 ), normal TE, 

\E\ N 

b p = —e 3 j— j v ' S[u] T = -e 3 0V ■ S[u] T , light weighted TE, 

I I fe= 1 

where $ = |£|-/V/|V| denotes the volume fraction. We use here the asymptotic approximation 
dvo/dt = du/dt = Re -1 V ■ S[u] T + Fr~ 2 e s + CJ(e 2 ) (cf. Theorem fill) and Assumption [2l (A8) 
that allows us to evaluate the gradient of the velocity at x instead of the center of mass of the 
corresponding particle. 

Remark 12 (Macroscopic definition of the local quantities). In the derivation of b p we used the 
law of large numbers to replace the discrete arithmetic mean of N particles with the expectation E[.], 
which is in consistency with Assumption \10\ This transformation was especially done to achieve 
a description that is easily comparable with the results presented in literature, see e.g. mnang. 
However, it involves a small technical issue: While the ensemble average needs only the well-defined 
velocity fields of each particle m\ in V(x, i), the expression E[dvi/dt](x, t) presupposes the definition 
of an underlying random velocity field Vi : fi x R^ x W —> R 3 whose characteristic for the kt.h 
particle v 3 (t,w) = Vi(Cq (t,u>),t,w) solves almost surely the corresponding ODE (14.Gal) . Having the 
meaning of the averages of the particle related quantities in mind we stick to the presented notation 
for reasons of readability. 

Prolate ellipsoids. For a prolate ellipsoid the lengths of the semi axes satisfy d\ > d 2 = d 3 > 0. It is 
convenient to introduce the aspect ratio a r = di/d 2 > 1 and the parameter v = (a 2 — l)(a 2 + 1) _1 . 
Its rotational behavior can be expressed in terms of the main director p 1; (Ro)ij = e; ■ p ; as the 
angular velocity w 0 (I4.6b[) becomes 


w 0 = 0.5V X x u 0 - v(p 1 ■ E[u 0 ] • p 3 )p 2 + ^(p 2 ■ E[u 0 ] • p 3 )p 3 
= 0.5V X x u 0 + nPi x E[u 0 ] ■ p 1; 


(4.9a) 


due to the orthonormality of {p 3 , p 2 , P 3 }. In combination with 


|pi=^oX Pl , 


(4.9b) 


m is often called Jeffery’s equation in remembrance of [10] . In the asymptotical context the 
equation was derived for normal tracer ellipsoids in the work of Junk & Illner El. As for the 
particle-stresses (14.81) . it can be shown by a technical but straight forward calculation [23] that 


±. r T 1 

-R 0 • M R 0 = a 1 p 1 ® p 3 O p 3 ® p x : E[u 0 ] + a 2 (p 3 <S> Pi • E[u 0 ] + E[u 0 ] • p 3 ® p 3 ) + a 3 E[u 0 ], 
0 


with d = did 2 and the geometry dependent M and cq as given in Abbreviation [I3l (cf. Appendix lA.il) . 
The expression holds for every particle k. Using Assumption [2] and shifting the isotropic part 
(M^ : diag(a°, a 2 , a 2 ))l into the pressure of the Newtonian stresses, we can state the particle- 
stresses in the well-known form with volume fraction <f> and deformation gradient tensor E[u] (see 

e-g- .13), 


E p = e 3 2(/)Y. P , 

i P = aiE[p 3 < 8 > p 3 < 8 > p x < 8 > p 3 ] : E[u] + a 2 (E[p 3 ® p 3 ] • E[u] + E[u] • E[p 3 ® pj) + a 3 E[u], (4.10) 
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Abbreviation 13 (Geometrical parameters of prolate ellipsoids for particle-stresses). The geomet¬ 
rical parameters a±, 02,03 of (14.101) are (cf. [ill 119]) 

4 


a 1 = 


61 = 


7i 


did 2 \4 72 0°d 2 


2 ~ 2&i — b 2 1 , 

2 

a 2 (2a 2 0 — 6—1) 
2did 2 (a)i — 1)02 (d 2 o 2 + ^i a i) 


— ^2, 


a 2 = 


bo — 


■ In 


a r + (a 2 — l) 1 / 2 
2a r {a 2 r - 1)V2 ^a r -(o2-1)1/2’ 
1 


X = 


46i 

did 2 ’ 

1 

Wf’ 

2a r 6 


03 = 


46 2 

dldn 


did? a2 — 1 


(a 2 0 — l) 


1 




1 


didf 4(a2 - l) 2 


71 did2 4( a 2-l) 2 
In addition, 


(30 + 2 a 2 - 5) , 

(—(4a 2 — 1)9 + 2a 2 + l) , 


did 2 aj! — 1 

1 1 


7 2 = 


did| (a2 — l) 2 
1 a 2 
did| (a2 — l) 2 


(-0 + 1 ), 

(-3 a 2 r 9 


2) 


((2a 2 + 1)0 - 3) . 


where the coefficients are given by 

A = — 

72(2^22 — ^33) — 7 i^n 



A 

H 

G* 

M = 

H* 

B 

F 


G 

F* 

C 

2 7 iM 

11 — 

72(^22 ^ 


6(72 + 27i°7 2 °) 


B = - 

F = - 
G = — 
H = - 


6( 72 o2 + 2 7l 0 72°) 

£* 2(^23 + A 32 )/2 — dl0<j{A 32 — A 23 )/2 

4/3° ^ 2^2 

<^ 2(^13 + A 3 i)/2 — d\(3 2(^13 ~ A 33 )/2 
2/3 2 (df a? + dla° 2 ) 

c*i(7li2 + A 2 i)/2 — dlPZ(A 21 — Ai 2 )/2 
2/3 2 (dfa? + d|a^) 


C = — 


F* = — 
G* = - 
H* = - 


72(2^33 ~ ^22) — 7 i-Aii 
6( 72 o2 + 2 7 l ° 7 2 °) 
^ 2 ( 7^23 + A 32 )/2 + d20i(A 32 — A 2 3)/2 
4^°d^a£ 

Q°(t1i3 + A 3 i)/2 + d 2 /? 2 °(Ai 3 — A3i)/2 
20 2 {dla1 + dlaf) 

02 ( 7^12 + 7l 2 i)/2 + df/3 2 (A 2 i — A V2 )/2 


2PZ(d\a 0 1 + d 2 2 af) 


(cf. Lemma\j£ 1/ with A zj = (A x ) ijt Ai = zaB((0, -p x ■ E[u 0 ] • p 3 ,Pi ■ E[u 0 ] • p 2 )) - Rq • E[u 0 ] • Ro- 

Corollary 14 (Asymptotical model of a suspension with prolate weakly inertial tracer ellipsoids). 
Let the assumptions of Theorem \ll\ be valid and let the volume fraction (f) = |£|iV/|V| be independent 
of x, then the suspension of prolate ellipsoidal particles is macroscopically described up to an error 
ofO(e 4 ) by 

Re (<9iU + (u • V x )u) = ? ? V X • S[u] T + e 3 2<^V • ± pT + f p + ReFr" 2 e g , V x • u = 0, infix R+, 

~ p 

supplemented with appropriate initial and boundary conditions. It is E according to (14.101) and 

ReE [^vi] , heavy TE, 
else. 


f 1 + e 2 (j) -pj — e 3 </>, 

heavy TE, 

( 

b + tV)®- 1 ), 

normal TE, 

f p = l 

(1 - "Y 

light weighted TE, 

l 
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Remark 15 (Properties of additional inertia related forces in suspension model). 

1) Compared to the classical results in literature, e.g. he our model includes a change in the 
overall Newtonian stress of the fluid and additionally, in the case of heavy tracer ellipsoids, 
an extra body force. The source of these effects can be identified in (I3.5cl) to originate from 
the one-particle associated functions ff for i = 1,2, which in turn have the role of source 
terms in the ODEs for c k according to the solvability conditions (|2.11l) . clearly seen in the 
case of ellipsoidal particles (14.6al) . This implies that the additional macroscopical effects 
have their origin in the deviation of the particle center of mass from the streamlines of the 
undisturbed fluid, which is given by c k — Cg ss ec k + e 2 ^". This dependence is also clearly 
observed for ellipsoidal particles. In the general case the solvability conditions for the linear 
and angular velocity may not decouple, but still the f, terms are associated to the ODE- 
system for the correction of particle movement. Thus throughout this text we refer to these 
terms in Theorem m as originating from the small disturbance. 

2) In the special case p = 1, we see by (12.41) that m/\£\ = 1 and the contribution due to the 
relative motion vanishes, similarly to the classical results in [T . 

Remark 16 (Splitting of the inner force). .4s mentioned in Section 1 3 . 11 the splitting approach 
(ED for the overall inner force into a force generated by surface stresses and a volume force is only 
needed, since the linearity of the mean value with respect to differentiation is not generally valid for 
volume-based averages. Here, we want to illustrate this fact based on our results: As we see from 
Theorem m the volume average for the surface stresses yields (apart from the Newtonian stresses) 
only Z p that can analytically be simplified to the classical particle-stress term of Corollary \lff This 
is a contribution of 0(e 3 ) which is dependent neither on the local coordinates y of a single particle 
-and thus cannot generate an effective force of 0(e 2 )-, nor on the mass function ot mass -and thus 
cannot change despite different inertia models. In contrast, the average of the volume force E[b] v 
generates a contribution of 0(e 2 ) and changes accordingly to the choice of the inertia model. 

5. Conclusion 

In this paper we presented a model for a suspension of small particles, which involves inertial 
effects. The latter were identified as the ability of particles to deviate from the streamlines of the 
surrounding fluid as a consequence of different scaling of the particle momentum balance. This 
was achieved by an appropriate scaling of the particle mass, respectively the density ratio. To 
keep the results of HE we retained the basic asymptotic approach and restricted our choice of 
the mass function accordingly. Having characterized the microscale behavior of the particles, we 
modeled the particle suspension following H] , supplementing a strategy to take additional forces into 
account, which have their origin in the relative motion of the particles. Afterwards we gave different 
models for the corresponding inertial particle regimes. These models are composed of incompressible 
Navier-Stokes-like equations with modified, non-Newtonian stress and force. Besides the classical 
part in the stress, we found a modification entering the overall viscosity of the fluid. To illustrate the 
general results, we applied the theory to the classical example of symmetrical ellipsoidal particles. 

One practical restriction of our approach is the fact that we cannot give a recipe on how to 
choose an inertial regime for a given particle and fluid, since our model of inertia was formulated 
by means of an asymptotic behavior of the particle. For a concrete situation, one has to observe 
the particle motion and then, based on the intensity of deviation of the particle motion from the 
fluid streamlines, decide which regime is applicable. 

Appendix A. Analytical statements for ellipsoidal geometry 

In the appendix we summarize some fundamental results of different works for ellipsoidal particles, 
the individual results are from Oberbeck HE Edwardes [B], Jeffery HHj and Junk & Illner HD- Since 
we use especially the solutions of Oberbeck and Jeffery in the derivation of the suspension model 
for ellipsoidal particles in Section [4j we provide here the relevant arguments. In Appendix IA.1I 
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we introduce all relevant functions needed to formulate the solutions of Oberbeck and Jeffery in 
Lemma IA.21 also analyzing the decay properties of the involved quantities in Lemma IA.1I In 
Appendix IA.2I we briefly present the necessary steps for the analytical computation of the surface 
moments arising in the solvability conditions of Lemma [5] by using Lemma IA.2I 


A.l. Oberbeck and Jeffery solutions. As stated in Section [2 an ellipsoid is a set £ = DBi 
with D = diag(di,d 2 ,cfe), di > 0 and unit ball B\ in R 3 . This implies £ = {y € R 3 | |JD -1 y|| 2 < 1}. 
Consider the function A : R 3 \ £ —> Rj defined as 


d\ + A(y) + d 2 2 + A(y) + dj + A(y) 


(A.la) 


whose existence and regularity are guaranteed by the implicit function theorem. The formulation of 
solutions for the Stokes problems in R 3 \£ are based on the following geometry associated functions 


<5(A) = det(D A ) 


pOO 

1/2 , XW=J <5(s) _1 ds, 

poo 

Pj{ A) = / (d 2 j + s)6(s)~ 3 ds, 



(d 2 + s) 1 <5(s) 1 ds, 
(d 2 + s)s5(s)~ 3 ds 


(A.lb) 


for j = 1,2,3, with D A = D a (A) = diag(d 2 + A,* = 1,2,3). When a 7 -, (5j , 7 j and y are eval¬ 
uated at zero, we abbreviate the function value with the index °, e.g. y° = J 0 °° <5(s) _1 ds. Ad¬ 
ditionally, we introduce ipj( y) = /3y(A(y))j/fc, where ( j,k,t ) is a permutation of (1,2,3), and 
w(y) = S(s)~ 1 f e i( y, s) ds, where f e i{ y, A) = y ■ D)) 1 ■ y — 1. In the following some derivatives of 

the above functions are needed. We set n( y, A) = (y ■ Di) 2 • y) _1 , in other words d\f e i = —(i _1 . It 
follows for y g R 3 \ £ 

5yX'(A(y)) =-2^y ■ D” 1 , dyui = 2y ■diag(a 1 ,a 2 ,a 3 ), (A.2a) 

9 yy x(A(y)) = -2^ (D^ 1 - 2,i (D^ 2 ■ y ® Di) 1 • y + Di) 1 ■ y ® Di) 2 ■ y)) (A.2b) 

- 2^ (- tr(Di) 1 ) + 4/iy • D^ 3 ■ y)D" 1 • y ® Di) 1 ■ y) , 
d yy w = 2diag(ai,a 2 ,Q!3) -4^ (D)) 1 ■ y <g> D)) 1 ■ y), (A.2c) 


, /I f djkUi $ijVk “b dfk!J'j j 

^ \(^f + A )(d 2 + A) (d 2 + A)(d 2 + A) J 


(A.2d) 


0 /j 2 ViUjUk {{d 2 + A) 1 + (dfc + A) 1 + (d 2 + A) 1 + 0.5tr(D A 1 ) - 2/ry ■ P A 3 ■ y) 
5 (d 2 + A)(d 2 + A)(d| + A) 


Let (i,j, k ) be an even permutation of (1, 2, 3), the derivatives of %/ii read as: 


— 2 jg (D A )^(D^ ' Y ^UjUk "b PiidljUk “b (A.2e) 

dymvi^i = ~b dikdjm) ^UjUk^ ^ (2/r(D^ ■ y) m (D^ ■ y)i J- (D A )j2(D^ )^m) (A.2f) 

- 2^(D A )jj ((D)) 1 ■ y )m(Sejyk + SekUj) + (Di) 1 ■ y )e(Sj m yk + SkmVj )) 

u 2 

+ 4»ji/fc^-(D A ) ij ((D)) 2 ■ y^D^) 1 ■ y) m + (Di) 2 ■ y) m (D^ 1 ■ y) t ) 

— 21/jT/fc —(D A )ii(4/ry ■ D A ‘ - y — 3tr(D A ))(D A • y)^(D A • y)m- 
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From (IA.2bl) . (IA.2cl) and (IA.2fl) it follows 


AyX 

Ay ifi 


Ay LO 


-2 J(tr(D* *) - Apy ■ D^ 3 • y + (— MD/ 1 ) + Apy ■ D/ 3 ■ y)) = 0, 

~4—(D^)jjj/y2/fe(tr(D^ ) — (D a )»») — 2-^yjyj-(2 + (D A )ji tr(D_^ )) 

2 

+ 8^(D x ) iiyj y k y ■ D/ 3 • y - 2^(D A ) ii t/ j y fc (4/zy • D/ 3 • y - 3tr(D^ 1 )) 


4 -2. 


tr(D: 3 ) 


A(y) 


d » + J =- 4 


'A(y) 


a,j d * + j)=o. 


(A.3a) 
(A.3b) 


0, 


(A.3c) 


since d\5 1 = —0.5<5 1 ti/D,^ 1 ) and liin^oo S 1 = 0. This is an important feature of the functions, 
which will be used in Lemma lA.21 Next, we present some asymptotic properties of these functions 
for ||y|| — > oo. 

Lemma A.l (Asymptotical properties). The geometry associated functions defined in (IA.1I) fulfill 
the following properties for r = ||y||, r —> oo: 



where the notation f ~ </> stands for f = <f> + o(4>). 
Proof. From the definition of A(y), it holds 


i = y ■ d a 1 • y 


< (mini d 2 + A) V 2 
> (max* d 2 + A) _1 r 2 


and thus 0 < minj d 2 < r 2 — A < max.i d 2 , which implies 

, r 2 - A 

lim --— = 0. 

r—>■ oo r A 


Set d = max.i di , then the behavior of the algebraic functions S and p results as 



For the integral functions the following statement is used: Let g : (0, oo) —> R + be a continuous 
function with gds < oo. Define f(t, Q) = f® gds, then for any Q £ R + it holds 


/(A(y), Q) - f{r 2 {y),Q) 


r r2 ( y) 

/ gds < (r 2 
J A(y) 


A) max q < d 2 max q. 

se[A,r=] se[A,r-2] 


Since Q was arbitrary, this is still true for Q —> oo. Another statement also needed is: Let a, a, s > 0 
then 1 — (a/s + l) -a < aa/s. This can be directly concluded from the fact that the function 
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/(s) = aa/s + ( a/s + 1) “is strictly decreasing and lim s _ i>00 / = 1. Consequently, 


X- 2 r 


-l 


2 r 


l 


1 

— 2r _1 

d 2 


S(s) 1 ds 


poo 

/ (5(s) _1 — s -3 / 2 ds 

Jr 2 

d 2 i r°° 

< -- max (i(s)^ 1 H -- / s~ 3 ^ 2 — (5(s) _1 ds 

“ 2r _1 s e[A,r- 2 ] v > 2 1 — 1 J r 2 v ’ 


d 2 \~ 3 / 2 l 

<-— + 


.-i 


2 r 

d 2 2 
< — r“ 2 
“ 2 


2r _1 

-3/2 


W 3 / 2 ( 1 - 


(A-‘ +1) 


-3/2 


ds 




0. 


The same steps lead to 


a»— 2/3r 


-3 


2/3r - 3 
A - 2/5r -5 


2/5 r - 5 
7t - 2/3r -3 


2/3) 


.-3 


< ——r 

“ 2/3 


< 

“ 2/5 

d 2 2 

< -^r“ 2 
“ 2/3 


-5/2 

“2 ' + 1 


-7/2 


+ 1 


-5/2 


0. 


□ 

We summarize some results of DU El ESI for the disturbance flow of an ellipsoid in a Stokes flow 
in Lemma TA. 2 1 

Lemma A. 2 (Oberbeck and Jeffery solutions). Consider an ellipsoid with its geometry associated 
functions. Let v £ 1R, 3 . A £ ]R 3x3 be a constant vector, respectively matrix, where tr(A) = 0. Define 

uob = (x 1 - VyX ® y + O.bdyyU} • D 2 ) • O v, (A.4) 

(Uif i\ 

(A.5) 



d Vl ifi d yi tp2 d yi ^3 




/ t/Ai \ 

Uj e = 

dy 2 lfl A /2 ^2 dy 2 1p 3 



H - Vy X 

VV>2 


dy 3 l/j 1 A,3^2 A/s'As, 


w 


\Wlf3J 


dyyCO M t y — M V 


yUJ, 


as well as pob = —2V y x • 0 ■ v and pj e = 2 M : d yy u> with the matrices 
O = diag(x° F d 2 cx° ,i = 1, 2, 3) 1 , M 


A H G* 
H* B F 
G F* C 


and the corresponding coefficients given by 

A _ 27 /An - 7 /A 22 — 73^33 

6(7273 + 7i 73 + 7i7/) 

27 IA 22 — 73 A 33 — 7 i An „ _ 27 gA 33 — 7 /An — 7 /A 22 

‘ 6 ( 7 / 7 / + 7i°7 3 ° + 7i°7 2 °) ’ ~ ~ 6 ( 7 / 73 ° + 7i° 7 3 ° + 7i°7/) ’ 

p, _ Q 2 (^23 + A 32 )/2 - rf§/3/(A 32 — A 23 )/2 _ ct/(A 23 + A 32 )/2 + d|/3/(A 32 - A 2 3)/2 

2B?(d%a% + d?a?) ’ 2/3/(do a/ + d? a/) 


B = 


„ _ _ ^ 3(^13 + A 3 i)/2 — d\j 3 2 (Ai 3 — A 3 i)/2 _ _ a:/(Ai 3 + A 3 1)/2 + d 3 /3/(Ai 3 — A 3 1)/2 

2 /3/( d 3 a 3 + d l a l) ' 

A 2 1)/2 — d|/3 3 (A 2 i 
2/3/(ci?a/ + d|a§) 


2 /3/(d|a/ + dfa/) ’ ^ 2 /3/(d|ag + d?aj) 

a/(Ai 2 + A 2 i)/2 — d|/3 3 (A 2 i — Ai 2 )/2 * a/(A 32 + A 2 i)/2 + d 2 /3 3 (A 2 i — A 32 )/2 

11 =-—;—r-^-— ;-, ii =--— ---:- 


2/3/(rf2 a o +c ?2 a o) 
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R — (^23 + -A 32 )/(2/3°), 

U = 2d\B - 2d\C, 


S = (A 13 + A 31 )/(2P° 2 ), 
V = 2djC - 2 d\A, 


T = (A 12 +A 21 )/(2p° 3 ), 
W = 2d\A - 2d 2 2 B. 


Then (u,p) G {( uob,POb ), (uj e ,Pj e )} is a solution of 

Vy ■ S[u] T = 0, Vy ■ U = 0, 

u = h, 


y G R 3 \ E, 

y e d£, 


u — > 0, ||y|| -MX), 

with the decay properties ||u|| < c||y|| 1 and ||V y u||, \p\ < c||y||“ 2 , c, c > 0. The Dirichlet condition 
reads as h = v for u = uob and h = A • y for u = uj e . 


Proof. The proof of this lemma can be found in [HI for the translational boundary condition and in 
pTO] for the linear one. For the sake of completeness, we show the relevant steps here. First, using 
0 = AyX = AyW = Ay if>j from (1A.3I) implies 

Vy • UOb = (VyX - VyX - ^ Ay* + 0.5VyAyW • D 2 ) ' 0 ‘ V = 0, 

Vy ■ Uj e = ^ Ay Ipjaf + VyAyW • ' y + M : DyyOJ — M : dyyOJ = 0. 

j 

by means of the identity V • V x w = 0 for any vector w. Here, a R = ( R , S, T) T . Moreover, 

AyUOb = (AyXl — VyAyX ® ~ 2Dyy\ A 0.5<9yyAyU 1 ■ O V = Vy( —2Vy% O V), 

Ay U,/ e = Vy ^2 Aylpjaf + Vy X (U Ay ^1 , V Ay lf 2 , W Ayll) 3 ) T + 8yy AyU> ‘ y 

j 

+ 2Vy9yyW ) — M ' VyAyW = Vy(29yy UJ \ M). 

The matrices O, A and the coefficients R,...,W are chosen in a way such that that the boundary 
conditions on 8£ hold: 


D 2 j O v 

= (x°l + diag(a°, i = 1,2,3)- D 2 ) • diag(x° + d 2 a°, i = 1,2, 3)" 1 v = v, 

analogously for uj e . Last, the decreasing properties for r = ||y|| —> oo are shown. This follows from 
the definition of the velocities (IA.4I) and (IA.5() as well as the corresponding pressures, 

||u 0 b|| < c(|x| + ||5 y x||r + ||9yyw||), \p 0b \ < c||9 y xll, 

||uj e || < c(max||dyV’fc|| + ||<9yyw||r + ||d y w||), \p Je \ < c||<9y y w||, 

k 

||9yUOb|| < c(||5 y x|| + ||<9yyX|k+ ||9yyyW||), 

||9yU Je || < c(max||dyy^fc|| + ||<9y yy w||r + ||<9 yy u;||), 

k 

with appropriate constants c > 0. For r sufficiently big, it is known from (IA.2I) . Lemma [A. 1 1 that 
Ixl < cr" 1 , ||6»yXll < cr~ 2 , ||6» yy x|| < cr" 3 , ||<9 y w|| < cr ~ 2 

||<9 yy w|| < cr -3 , 11 dyyyOJ 11 < cr" 4 , \\dytp k \\ < cr" 4 , \\dyyif k \\ < cr~ 5 , 

yielding the proposed decay properties: 

||u 0 b|| < cr" 1 , |po&| < cr" 2 , ||<9 y uob|| < cr" 2 , 

||u. 7 e|| < cr" 2 , |pj e | < cr" 3 , ||i9yUje|| < cr" 3 . 

□ 


u Ob — X°l + 2 


(y-D -y )- 1 

det(D) 


D • y ® y, 


+ 0.5 ( 2diag(a°,i = 1,2,3) - . y D" 2 • y 
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A.2. Determination of surface moments. For the solvability conditions (12.111) we provide the 
surface moments of an ellipsoid. Computing the Newtonian stresses of Oberbeck and Jeffery solu¬ 
tions for the Stokes problems fLemma lA.2l) on d£ yields 

S[uob](y) = y (-(v ■ D -2 • y)l - v ® D~ 2 • y - D 2 ■ y 0 v + 2 ■ D 2 • y)D 2 ■ y 0 D 2 • y) , 

,,2 

S[uj e ](y) = A + A T — 16^r(y D 2 • M D 2 • y)D 2 • y 0 D" 2 • y 

0 

+ “F (M • D 2 ■ y 0 D 2 ■ y + D~ 2 • y 0 M ■ D 2 • y) 

0 

- 4M : (diagK, a° 3 ) - 2^D- 2 • y 0 D" 2 • y)l, 

where p- = { y ■ D 4 ■ y) _1 , (5 = det(D) and v = O ■ v. For z € ffi 2 the following identities hold: 

S[uob](D ■ z) • D 1 ■ z = —y, (A.6a) 

0 

S[u Je ](D ■ z) • D 1 ■ z = ( A + A t + - 4(M : diag(a°,* = 1,2, 3))l j ■ D' 1 • z. 

The last result simplifies further, if A is skew-symmetric, i.e. A = f?(v) for some v £ It 3 . Then, 

S[u Je ](D ■ z) • D 1 ■ z = -^S(diag(ci,i = 1,2,3) • v) • D • z, (A.6b) 

with the constants d = (V 2 a° + V 2 a£) _1 , where (■ i,j , k) is a permutation of (1, 2, 3). 

To evaluate the integrals over the surface of an ellipsoid we transform them to the unit sphere. 
According to tni it holds: 

/ F(y) • n(y)ds(y) = f F(y) ■ °_ 2 V - da(y) = [ F ( D ■ z) ■ °_i * det(D)||D~ 1 ■ z||da(z) 
Jas Joe || D ■ y || V$? || D - z|| 

= 6 I F(D ■ z) • D 1 • zds(z), (A.7) 

Vs? 

for an integrable function F : R 3 —>■ R 3 X R 3 . 

Now consider the six Stokes problems in Lemma [5] According to Lemma IA.2I the analytical 
solutions (w q,p q ) are given by (IA.4I) with v = e g for q = 1,2,3 and by (IA.5I) with A = B(—e q - 3 ) for 
q = 4,5,6. To calculate the surface moments s q ,t q , V g and W q , we first apply the transformation 
to the unit sphere (IA.7I) and afterwards use (IA. 6 al) or (IA. 6 bl) for the corresponding moments: 


s q= [ B (y) ■ S[w g ] • n ds(y) = -4 

[ V?(D ■ z) • Vq(z)ds(z), 

Jd£ 

/ s? 

t q = f S[w g ] • nds(y) = -4 

/ v -j( z ) ds ( z )> 

V d£ 

/ \ 

( v ?)ii = ( [ 2/*S[w g ] • nds(y) N ) = -4 

( / (D ■ z);v g (z)ds(z) ) 

\Jd£ Jj 

\M Jj 


(W q )ijk = ( [ yiVjS[w q ] • nds(y)) = -4 (I (D z) f (D ■ z) 3 ■ v 9 (z)ds(z) 

\Jd£ J k \Vs? 

with M q = O e q and v 9+3 = —c q B( e q ) ■ D • z for q = 1,2, 3. Since the integral of a homogeneous 
polynomial p{ z) fulfills f s 2 p(z)ds(z) = 0 if p has odd degree, it follows 

s q = 0 , M q = 0 , for < 7 = 1 , 2 ,3, 

t q = 0, W q = 0, for q = 4, 5, 6. 
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The other moments follow, considering / S 2 (D • z)$(D • z)^ds = (D 2 )^-47t/3 and |S 2 | = 47r, 
(Sqf+3)fc = ^ ^ ^ )im = ^ ^ ^ (fiq&q)z {$kZ$im $km$iz){^ 


i,j,t,m 


iim 


= -c^. ((D 2 - tr(D 2 )l) e q ) k = C 9+3 |£| ((D 2 - tr(D 2 )l) e q ) k , 

i^q+s)ij = g ^ )it e kej{ c q e q)k = ~ c q g (D • B(e q ))ij = Cg+ 3 |£|(D • B(e q ))ij, 

k,l 

t q = -167tO ■ e q = 3£ g |£ \e q , 

(W q )ijk = D,y (O • e q )k = Cq|^|(D <8 e q )ijk, 

for q= 1,2, 3, where £ q = —4det(D )~ 1 O qq and £ g+ 3 = —4det(D) _1 c 9 , with O given in Lemma liOl 
and the constants c, = {fa° + d^c^) -1 , here (■ i,j, k ) is a permutation of (1, 2, 3). 
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